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We investigate the non-linear evolution of spherical density and velocity perturbations of dark matter and dark 
energy in the expanding Universe. For that we have used the conservation and Einstein equations to describe the 
evolution of gravitationally coupled inhomogeneities of dark matter, dark energy and radiation from linear stage in 
the early Universe to the non-linear one at the current epoch. The simple method of the numerical integration of 
the system of non-linear differential equations for evolution of the central part of halo is proposed. The results are 
presented for halo of cluster (fc = 2 Mpc~^) and supercluster scales (fc = 0.2 Mpc“^) and show that quintessential 
scalar field dark energy with small value of effective speed of sound Ca <0.1 can give noticeable impact on the 
formation of large scale structures in the expanding Universe. 
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INTRODUCTION 

Dark energy is mysterious dark component responsible for accelerated expansion of the Universe became 
an object of numerous studies in the last two decades. There are many possible explanations m, m of the 
nature of dark energy. One of the most promising among them is scalar field dark energy, which can be 
modelled as perfect fluid. In fact, having just a few parameters - the density p^e = Pcr^de) th® equation 
of state (EoS) parameter w^e = p/p and effective speed of sound = 6p/6p, this model can perfectly 
explain latest cosmological observational data m, ca), keeping a large volume of parameters space of 
possible values. Therefore, it would be interesting to study behaviour of perturbations of such dark energy 
on non-cosmological scales. This may yield more constrains on parameters of the model. 

It was already shown [2j, |15) that non-linear perturbations of scalar held dark energy can in principle 
influence structure formation even on galaxy scales. Finally, in the last decade there was a number of 
papers (see [3j for review), that studied accretion of dark energy on compact objects (black holes mainly) in 
hydrodynamical approach. In our works CD El we have shown that nevertheless such stationary accretion 
does not change the gravitating mass of the central object, in principle, it can give detectable contribution 
to dynamics of bodies in its vicinity, which would be possible to extract from current observations. 

Here we analyse the dynamics of dark energy in the collapsing halo of dark matter. In the next chapter 
we present the system of equations, which we use to solve the problem of evolution of spherical scalar 
perturbations in 3-component expanding Universe from early linear stage to highly non-linear one, when 
dark matter halo forms. In the last chapter we present numerical solutions and discuss their main features. 


EQUATIONS FOR EVOLUTION OF PERTURBATIONS 

We suppose that the Universe is spatially flat and filled by matter (m), dark energy (de) and radiation (r), 
the metric of the background space-time is Friedmann-Robertson-Walker (FRW) one. Each component is 
described in the perfect fluid approximation by energy density e, pressure p and four-velocity u*. The equation 
of state for each component can be presented as p = we, with Wm = 0 for matter, Wr = 1/3 for radiation 
and Wde < —1/3 for dark energy. We assume that dark energy is scalar field with Wde = const. The goal of 
the paper is analysis the evolution of spherical halo from the linear stage in the early radiation dominated 
epoch, through the turnaround point to highly non-linear stage, infall of matter before virialization, at dark 
energy dominated epoch. The local spherical perturbation distorts the FRW metric, so that it becomes 

ds^ = - a\t)e^^^^’^'>[dr^ + r\d0^ + sin^ 9dip% ( 1 ) 

where the metric functions u{t, r)and p{t, r) vanish at cosmological background. At the linear stage, when 
nit, r) ~ p{t, r) <C 1 the metric jT|) becomes the metric of conformal-Newtonian (longitudinal) gauge [1] in 
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spherical coordinates, since ~ 1 + 1 /, ~ 1 + /r. We rewrite the energy-momentum tensor of dark 

energy, = (cde Pde)'>J‘i {de)'^de ~ ^iPde^ ^ terms of proper 3-velocity of fluid v^e (measured in local 

frame), which has only radial component. The relation between components of 4-velocity and 3-velocity 

Vde of fluid are as follows 
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where v^e is in the units of speed of light. The non-zero components of energy-momentum tensor are 
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We decompose both density and pressure on background averaged and perturbed parts as Ede = ^dei^ + ^de), 
Pde = £de{wde + c^^de — 3a(l -|- Wde){c^s ~ '‘^de) f Vdedr), where density perturbation of each component 

is defined as 6{t,r) = {e{t,r) — £{t))/e{t). We study the model of dark energy, for which both Wde and 
are constant. The integral term in the pressure decomposition comes from non-adiabatic part of pressure 
perturbation of scalar field dark energy (details can be found in papers [5l |6l |T0]). The presence of this 
term makes our equation of state non-barotropic, as we work in frame different from proper dark energy 
one and hence, relation pde = WdePde holds only for averaged parts of this component. Taking into account 
this decomposition, and keeping terms with and only, we obtain the following energy-momentum 

tensor components of dark energy: 


Tq — £de{^ + dde) + £de(l + Wde + (1 + c‘l)dde)Vdei 


Tq ~ ^de 
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TI = -Ede {wde + cl5de “ 3a(l Wde){cl - Wde) J Vdedr^ - Ede (l + Wde + (1 + cl)5de) 
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Hereafter, one can get the corresponding equations and expressions for dark matter and radiation components 
from equations and expressions for dark energy just by putting Wde = = 0 for matter and Wde = = 1/3 

for radiation. 

To find the evolution of density and velocity perturbations we use two conservation equations, which have 
general covariant form 

rj^k ^ 1 ^ n 

^ dx^ 2 ^ 

(for i = 0 it is continuity equation, for f = 1 it is motion equation). Substituting our energy-momentum 
tensor one obtains: 
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The continuity equation for background density, which we also use, is Ede + 3|(1 + Wde)£de = 0- 
To hnd the metric functions we exploit the Einstein equations 


R] - ts‘R = ^ 
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If we construct the equation G{ — = x (T/ — Tg) then we will see that at the linear stages /U = —i^. At 

the non-linear stage, the right hand side of this equation equals zero at the center of perturbation too, which 
gives again reduction to only one potential. In this paper we analyse the dynamics of dark matter and dark 
energy in the central part of spherical overdensity only, therefore we accept /i = —u approximation, that 
gives us the possibility to use only one (q) Einstein equation for determination of one metric function u{t, r): 
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Here we used the notations Hat = -|- + e^), where N denotes the type of fluid and “0” marks 

the value at current epoch. In the flat three component universe Qm + 12^ + ^de = 1- Eqs. have 

non-relativistic limit (Appendix), which in case of dark matter (c^ = Wde = 0) coincide with well known 
classical hydrodynamic equations and Poisson equation accordingly. 

We are interested in the cluster and supercluster scales of perturbations, for which <C 1, Vm ^ 1 in 
the bulk of object, while 1 in their halos. So, the equations can be essentially reduced by neglecting 

of the terms like, Wm, I'Vde nnd higher order terms. The terms with vVv must be kept, since they are 
important at highly non-linear stage. The final reduced form of these equations, which we use in the code, 
can be found in our upcoming paper |in) . 

Therefore, we have seven Ist-order partial differential equations for seven unknown functions 

Smit,r), Vm{t,r), Sdeit,r), Vdeit,r), Sr{t,r), Vr{t,r), ( 8 ) 

which can be solved numerically for given initial conditions. 

At the early epoch the amplitudes of cosmological perturbations of space-time metric, densities and 
velocities are small and the eqs. (I5])-([7|) can be linearized for all components. Moreover, we can present each 
function of (t, r) as product of its amplitude, which depends on t only, and some function of radial coordinate 
r, which describes the initial profile of spherical perturbation, which can be expanded into series of some 
orthogonal functions, e.g. spherical ones in our case. In particular, we can present the perturbations of the 
metric, density and velocity of N-component as follows 
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In the analysis of the evolution of central part of spherical halo we can decompose r-function in the Taylor 
series and keep only leading terms: 
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where fkir) = sin kr/kr. It gives the possibility to reduce the system of seven partial differential equation 
for unknown functions Q to the system of seven ordinary differential equations for their amplitudes 
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which is presented in |10| . 

We set the adiabatic initial conditions in the following way. To find the relations between amplitudes 
at some ainit 'C 1 when the scales of interested here gravitationally bound systems were essentially larger 
than horizon scale, ainitk~^ 3> t we have used the linearization of system ©-([T]), which has analytical exact 
solution for one component Universe. Obviously, that for cluster and supercluster scales it was early radiation- 
dominated epoch, when S> Em ^dey so, the matter and dark energy can be treated as test components. 
The amplitude of metric function H is defined by density perturbations of the relativistic component. The 
non-singular solution of the corresponding equations has asymptotic values at ainit P*”®* = — C, = 

C, h®®®®* = C/[4:ainitH(ainit)], where C is some constant (see for details |8]). The solutions of equation for 
matter and dark energy as test components give the asymptotic values for superhorizon perturbations at 

= 3C/4, h^®* = C/[^ainitH(aimt)], = 3(1 + w)C/A, = C/[Aa^nitH{a^nit)]. These 

relations contains only single constant C, which value specihes the initial amplitudes of perturbations in all 
components. Below we put C = 2.6 • 10“^ for cluster scales k = 2 Mpc“^ and C = 6.5 • 10“® for supercluster 
scales k = 0.2 Mpc“^ at time when a = 10“^^. 


RESULTS AND DISCUSSION 

In this section we present the results of numerical integrations of the described system of equations. For 
that we have designed the Fortran-95 program based on open dverk.f package for ODEs, in which time 
depended functions are replaced by scale-factor depended, since a(t) is well dehned for given cosmology. 
The cosmological parameters are taken as follows = 4.2 • 10“®, = 0-7, Qm = 1 — Dr ~ ^de, Hq = 70 

km/s/Mpc. 

The results of numerical integration of system of equations for time evolution of amplitudes of dark energy 
and dark matter density and velocity perturbations and potential for different scales and parameters of dark 
energy are presented in hgs. 1 a)-f). In these hgures, the solid lines correspond for dark matter perturbations 
and dashed lines correspond for dark energy ones. The dotted lines at all panels show the predictions of 
linear theory. For such perturbations non-linearity becomes noticeable already at a ~ 0.1. 

One can note, that in all cases and for all sets of parameters analysed here dark energy is subdominant 
always: density of its perturbations is few orders lower then of dark matter’s ones. The reason of this is the 
large absolute value of pressure which keeps it from collapse. Large values of effective speed of sound (that 
is, speed of (also) pressure perturbation distribution) forces dark energy to oscillate after perturbation enters 
the horizon (hgs. la and Id ). Both velocity and density perturbation of dark energy are smaller during all 
time of evolution, and the difference is bigger for smaller scales, as can be seen comparing Fig. la, lb and 
Fig. Id, le. 

Another interesting thing is how behaviour of dark energy depends on its parameters. From the hgs. a) 
and c) one can see, that EoS parameter Wde slightly changes the character of evolution of inhomogeneities, 
the initial amplitude of density perturbation is ~ (1 -|- Wde)- At the same time, the effective speed of sound 
Cs actually dehnes how fast the perturbation will grow. Comparing hgs. a) and b), or hg. d), e) and f) 
we see, that the closer is to zero, the larger amplitude of perturbations is reached at the hnal stages 
of evolution. With decrease of effective speed of sound, the pressure gradient of dark energy decreases 
counteraction to gravity. The oscillation became smaller and even disappear at some really small values of Cg 
and the perturbations grow faster, time evolution of velocity and density become more and more similar to 
the matter’s one. This fact is in good concordance with character of dependence on effective speed of sound 
during accretion of dark energy on compact object, that we studied recently m- As can be seen from hg. 
f), when Cs = 0 the dark energy and dark matter has always the same velocity and the gap between their 
densities is constant and ~ (1 -|- Wde) which is velocity factor in the continuity equation for dark energy ©. 

These two features, together with the fact, that dark energy background density is smaller than other 
components’ in the past, imply, that dark energy could inhuence signihcantly on the process of dark matter 
hallo collapse only having specihc parameters - very small value of effective speed of sound. The current 
cosmological observational data do not exclude such models of dark energy f|13|L 

It is also interesting to point out, that predictions of the linear (doted lines) and non-linear theories at 
smaller scale for density and velocity become distinguishable approximately at the same time (scale factor) 
as for larger scale, however, for potential, linear and non-linear theories diverge only at late times for smaller 
scale. The reason for is, that initial perturbations where taken different for these scales, to have approximately 
the same density perturbation at a = 1. For the same initial amplitudes the perturbations of smaller scales 
would collapse first. 
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Fig. 1: Evolution of the amplitudes of density perturbations of matter Sm and dark energy Sde (top panels), velocity 
of matter Vm and dark energy Vde in the units of Hubble one (middle panels) and gravitational potential (bottom 
panels) at the central part of spherical halo of cluster scale which is collapsing now. In the top and middle panels solid 
lines corresponds to matter, dashed lines to dark energy and dotted lines at all panels show the prediction of linear 
theory. r 
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APPENDIX. Newtonian approximation of equations for evolution of 
DARK ENERGY PERTURBATIONS 


Conservation equation Tq.^ transforms to classic continuity equation for dark energy: 


^de + ~ Wde)^de + 0, ^ | (1 -|- C^)Sde'V + [l + U) -|- (1 -|- C^)(5] yv' + 


2v 


= 0 


Combination of conservation equations and T^.j^ transform to classic Euler equation for dark energy: 

~ ^ ^ ^ ^ + (^ + " 3w)(Irfe ^ Q 

a V r J a l + Wde + {'^ + cl)5de 2o a(l -h Wde + (1 + cl)5de) 

If we consider the last two equations for dark matter (c^ = Wde = 0), we will see that they coincide with well 
known hydrodynamic equation of collapse [12] (equation 9.17) 

The non-relativistic approximation of Einstein equation ([7|) gives the Poisson equation for metric function 
V in the coordinates of FRW frame 


The Newtonian gravitational potential <I> = z^/2. 
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